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/ ? a s 

The definition and applications of the motion of 
magnetic lines of force are reviewed and illustrated,, 

First, the concept of such a motion is introduced as an 
aid to describing the evolution of a known magnetic field* 

It is next shown that a known velocity field in a perfectly 
conducting fluid with embedded magnetic field may be regarded 
as describing the velocity of the magnetic field line sphere, The 
consequences of such a flow on the field are then derived, 
including the properties of flux and line preservation and 
Cauchy's general solution. The discussion is then extended to 
fluids having a finite high conductivity, with accent on the 
role of X-type neutral points in the magnetic field. The review 
concludes with a brief description of existing theories on the 
behavior of such neutral points. 



( 1) Introduction 


Every theory in the exact sciences is based on a math- 
ematical formalism which derives its results from some given 
fundamental relations. In addition, however, it is often possible 
to formulate statements of a qualitative nature or, at least, 
with some intuitive content, about the subject described by 
the theory. Such statements are very useful, for they give 
insight into the subject's behavior and allow the deduction of 
general results without excessive calculation* 

The statement "in highly conductive fluids, magnetic lines 
of force move with the material" belongs to this class. In 
complicated geophysical and astrophy sical situations this is 
often the first guide to the evolution of the system being in- 
vestigated, from which one may continue, say, by devising simple 
(e.g. two-dimensional) mathematical models. 

Unfortunately, the concept of the motion of field lines 
has at times caused misunderstanding among those not familiar 
with its precise meaning and formalistic background. The purpose 
of this article is to review this meaning and background and to 
describe the limitations of the concept and some of its past 
applications. It is hoped that by following the exposition the 
reader will gradually become familiar with the subject and acquire 
the "intuitive feeling" for it on which much of its application 
depends . 

A great part of the theory of motion of field lines parallels 
earlier work on vortex lines in ideal fluids. This follows from 
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the basic equation for the development of the vorticity Q 
in such fluids 

|2-Vx(yx£)=0 (1-1) 

which is formally identical with the equation governing the 
development of the magnetic field B in perfectly conducting 
fluids 

— — V x (v x B) = 0 ' (1-2) 

“ 

The latter equation is readily obtained by noting that the 
first order electric field E* in the frame of the moving fluid is 

E* = E + (v x B) (1-3) 

and that Ohm 1 s law then is written 

j* - a E* (1-4) 

Taking the curl of (l-4) and applying Maxwell's equations 
yields equation (l-2), in the limit as a-# 00 , j'V^ 0 * 

An early use of the concept of perfect conductors is due to 
Maxwell (1873; art. 654) who showed that a perfectly conducting 
closed shell shields its interior from magnetic variations arising 
outside it. Thus when the conductivity of any material tends to 
infinity, the magnetic flux existing in it becomes "frozen in" 
and ceases to be influenced by external fields. Such ideas of 
"frozen flux" were expressed by early theories of superconductivity 
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which, however, did not consider a fluid medium and which 
were abandoned after the discovery in 1933 of the Meissner 
effect. The viewpoint that in astrophysics ionized gases can 
often be regarded as highly conducting fluids in which M the 
liquid is 'fastened' to the lines of force" was first advanced 
by Alfven (1942-b) who also used the expression that in this 
case the lines of force are "frozen-in" into the fluid (Alfven, 
1942-a) . The extension of the theorem of conservation of 
vorticity, derived by Helmholtz in 1858 (and 'independently , by 
Dirichlet) to the case .of hydromagnetic flow, probably occured 
at about that time, for Walen in his work on solar magnetism 
(1946, eq. 12) refers to it as "well known". The fundamental 
vorticity equation of Cauchy (l8l5) was similarly translated 
into hydromagnetic terms by Lundquist and Cowling (Lundquist, 

1951 , 1952) and many of the remaining details of the analogy 
were worked out be Maiersberger (1964). In the present work, 
hydromagnetic equations which originated in the study of vorticity 
will be referred to by their original names. 

The theory described here deals with highly conducting 
fluids; most of its applications, however, are to plasmas. In 
order to justify the extension of (l-2) to plasmas, two extreme 
situations will be considered. If the plasma is dense and 
dominated by collisions, fluid equations may be derived. Ohm's 
law is then considerably more involved than j(l-4) (Spitzer, 1962 
eq, 2-2l) , but with appropriate approximations only the pressure 
VPi (subscript i for ions) and the gravitational field V0 have to 



be taken into account (Spitzer, 1962). 


m 

i = *(,§* -A Vp -zi 7 ®) (u5) 

e 

If the electron density n^ is approximately constant or if p^ 
can be regarded as a function of it alone, the additional terms 
disappear when the curl is taken and (l-2) is obtained as before. 
In strong magnetic fields a may be anisotropic, but as long as 
the current flows M easily M (in the sense of eq. 4“3) in any 
given direction, the high-conductivity approximation can be 
expected to hold. 

At the other extreme one finds the "guiding center fluid" 
in which collisions are extremely rare. It may be shown, 
however (Parker, 1957a; Thompson, 1962, ch. 8, §7) that if the 
magnetic field is sufficiently strong to maintain the guiding 
center approximation, the electric field E* of eq. (l-3) vanishes 
in this case to lowest order in the ratio between the gyration 
time and the time scale of the macroscopic motions. The basic 
equation (1-2) then follows, so that the guiding center fluid may 
be regarded as a nearly perfect conductor in this connection . It 
should be added, also, that the slippage of lines of force in a 
dense plasma due to finite conductivity is much more readily 
analyzed than that occurring near the fringes of the magnetic 
field in a guiding-center plasma, in a region where the guiding 
center approximation becomes inaccurate. 

The motions discussed here will be relatively slow ones. 
Nonrelativistic transformations (e.g. equation 1-3) will be used. 
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the displacement current will always be neglected and wave 
phenomena will not be included. 

(2) Euler Potentials 

The magnetic field is solenoidal and therefore can be 
represented by two scalars. This representation may be achieved 
in several ways; for studying magnetic lines of force it is 
advantageous to introduce Euler potentials (Euler, 1770, §26, 
§49; Truesdel 1, 1954, §13) usually denoted by a and 0, with the 
property 


B rr Va x V0 = V x (a?s) (2-1) 

By this definition B is tangent to each of the families of surfaces 
a - con stan t and j§ = con s tan t . The line of intersection of any 
two such surfaces therefore defines a line of force characterized 
by the associated values of a and 0. Other commonly used methods 
of representing B afford no such direct way of evaluating lines 
of force. 

The use of Euler potentials also facilitates the introduction 
of tubes of flux. Consider a flux tube with a rhomboidal cross 
section, defined by the 4 lines of force (a, 0), (a + da, 0), 

(a? 3 + d0) and (a 4- da, 0 4 d0). The flux enclosed by the tube, 
evaluated in any plane normal to B, is 

d§ ^ J B* d_s - a V0 * d Ju - <j) ad0 (2-2) 

On two sides of the rhomboid 0 is constant and no contribution 
to the integral arises, while the two remaining ones contribute 
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(a + da)dg and -(adg). Thus 


d$ = dadg (2-3) 

If a and g are supplemented by a third spatial coordinate y, 
the flux through an area S on a surface y = constant is then 

§ = J* dadg (2-5) 

s 

In static fields, the magnetostatic potential may be 
chosen as y. It should be noted however, that in general it 
is not possible to stipulate that y varies only along the line 
of force , i .e . that Vy is normal to Va and vp. Such a stipulation 
implies B = e?Y (e some scalar function of position) , from which 
B • (V x B) = 0, a condition which is not in general met. 

The Euler potentials a and g each satisfy the first order 
partial differential equation 

B • 7f a 0 (2-6) 

and the possibility of expressing B by them may be proved (e.g. 
Phillips, 1933; §20) by using the existence proofs for solutions 
of this equation. The potentials are not uniquely defined, 
however: any independent pair of solutions to (2-6) leads to 

some choice of a and g. Given a pair of such potentials, an 
alternate choice (e.g. Ray, 1963) is h(a, g) and g(a> g) » provided 

(2-7) 


d(g> h ) = 1 
d(a,3 ) 
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Euler potentials have two major disadvantages. First of 
all, the representation is not linear and therefore the super- 
position of Euler potentials due to several sources is not in 
general valid. As a result, it is not possible to derive the 
functions <x and g analytically, except in some simple field 
configurations. One such case is provided by axisymmetric poloidal 
fields, such as that of a magnetic dipole. The vector potential 
A may then be taken in the azimuthal direction, leading to the 
choice 


0 = 0 


A ^ <x(r, 9) . 

— r sin 0 ~ 9 


(2-8) 


In fluid dynamics eq. (2-8) is often used to describe an 
axisymmetric fluid velocity field, with (-a) called the Stokes 
stream function (Mi lne -Thomson i960, Chapt. 15; O'Brien. 1963)* 
Superposition of potentials exists in this case, all of them 
sharing the same g. Another case for which a can be derived 
analytically (and superposed) is provided by purely toroidal fields 


B - V x r | = Vi|t x t 2 

Still another is that of a two-dimensional field, where JB 
neither depends on z nor has a component in its direction. One 
may then choose 

B = V x Ai = 7 a (x,y) x VZ 
z 
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The second and more serious disadvantage of Euler potentials 
is that in most magnetic fields lines of force are not closed and 
therefore attain infinite length in a bounded region of space. 

Consider (McDonald, 1954) a current-bearing circle r - R in 
the z = 0 plane of a cylindrical coordinate system. Because of 
axial symmetry, the a surfaces are toroidal rings and the g sur- 
faces meridional planes in which the lines of force, which are 
closedjlie. Now let a current filament be added along the z axis, 
giving the lines of force a slant and thus transforming them into 
helices. It may be shown that the a derived without the current 
filament may still be used; however, a g surface now resembles 
an archimedean screw with its axis twisted around the circle and 
its ends, in general, not meeting. Except for those cases in 
which a line of force returns into itself after a finite length 
(there is an infinite number of such cases, but of measure zero), 
a line of force will range all over its ot-surface, with different 
sections coming arbitrarily close to each other but never meeting. 

Given a point on the cr“Surface it is possible to construct a 
g-surface through it. This surface, however, if continued for a 
sufficient length, will fill a region in space and approach 
arbitrarily close to any given point in it. It is then impossible 
to construct a f ami ly of g -surfaces continuously dependent on a 
parameter. The characterization of lines of force by families of 
Euler potentials is therefore not feasible when such lines are 
bounded but of infinite length . 

The difficulty may be lifted by introducing a surface of 
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discontinuity which artificially terminates the (a? B) lines 
(Dungey, 1958; §3°l)° In the example given the part of the 
plane z = 0 with r£R may be chosen as such a surface , with 
each pair of parameters (a» B) defining a segment of a line of 
force from the point where it leaves the top of the plane to 
where it reaches the bottom 0 One may still consider infinitely 
long lines of force but such a line is then characterized by an 
infinite number of p values, with (3 changing every time the line 
crosses the surface of di scontinui ty c 
( 3 ) Kinematics of Lines of Force 

The first problem associated with motion of lines of force 
may be posed as follows 0 Given a magnetic field which changes 
in time in a known way, can an alternative description of this 
time ^dependence be made by assigning a velocity to its lines of 
force? 

In classical fluid dynamics, the rate of change in some 
quantity X, as observed by a particle moving with the fluid's 
velocity is its '%iaterial derivative” 


dX _ 3 X 
dt = St + 


u ° 7X 


(3-l) 


If X is convected with the flow, the material derivative 
vanisheso It therefore seems reasonable to call u the velocity 
of the lines of force if it satisfies 


If + Si • - 0 


+ u • VB = 0 
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(3-2) 



A velocity satisfying the preceding condition can always 
be found (Northrop, 1963) > namely 

-t? ve ) x ? (3 - j) 

as may be verified by substitu tion . This solution is not 
unique: any velocity parallel to B may be added to it without 
changing the validity of (3-2), since B is normal to both Va and 
ve. Because the component of u parallel to B is arbitrary, it 
i s customary to set it equal to zero ; when di sen s sing the velocity 
of a line of force, then, only the component normal to B is 
implied . 

This definition (3“2) for the velocity of the lines of force 
is not entirely satisfactory: the Euler potentials are not 

directly observable and it is not a-priori certain that equations 
(3-2) have to be satisfied. Indeed, since at any time there exist 
many equivalent choices of a and g describing the same .B, one may 
always add to the velocity satisfying (3-2) another velocity 
which describes a time-dependent transformation of Euler potentials, 
without altering the physical picture. Such a velocity will be 
termed a "relabeling velocity". For instance, in an axisymmetric 
field one can add any (not necessarily steady) angular velocity 
around the symmetry axis, for such a velocity merely relabels the 
lines and has no observable effect. 

Following Newcomb (1958) we therefore broaden the definition 
and regard as a "velocity of lines of force" any velocity v 
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satisfying eq, (l-2) 



“ ~ V X (v x B) =0 
a t 

It may be verified that u of eq 0 (3-3) is one of the solutions. 
Newcomb showed that by this definition all "verifiable con- 
sequences" which are usually associated with the concept of the 
motion of lines of force are satisfied; the two most important 
of these consequences 3 namely line preservation and flux pre- 
servation 3 will be discussed separately in Section 4« 

The most general solution of (l-2) will be the sum of the 
particular solution (3-3) and the general solution w of 

V x (w x B) - 0 (3-4) 


If one specifies* as before* (w ° B) = 0 , this general solution 
has the form 


w 


x B 


(3-5) 


where tjj (a s 3) is a scalar conserved along a line of force „ 

Since for any "relabeling velocity" which expresses no variation 
in the field* equation (l 2) reduces to (3-4) & it is clear that 
such a velocity will have the form (3-5) • 

(4) The Motion of Perfectly Conducting Fluids 

In a conducting fluid the presence of an electric field 
causes current to flow. If the fluid is moving* the relevant 
electric field is the one taken in a frame of reference moving 
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with the fluid; to first order this field is given by eq. (l-3) 


E* = E + (v x B) 

In the simplest case the relation between E* and the current 
density j* in the moving frame follows Ohm's law and is linear 
and isotropic as in (l-4) 

j* = cr E* 

For ionized gases in a magnetic field the relation may be 
anisotropic (the conductivity q is a tensor) or even nonlinear, 
but in general V' and E* increase together. For simplicity, 
a simple ohmic relationship with scalar will be assumed here, 
as well as those approximations valid when v < c. By Maxwell's 
equation s 


—V X (v X B) = V X— (V X B) (4-1) 

St - | x Q a 

or, if <j is constant everywhere 

•sD 

-f-Vx(vxB)= — P 2 B (4-2) 

The case of interest here occurs when o is very large, so that 
the term involving it may be neglected and (4*l) reduces to 
( 1 -2) . A rough but useful criterion to indicate when such an 
omission is permissible is furnished by the magnetic Reynolds 
number (Elsasser 1955> 1956, 1957). If L is a typical dimension 
over which B varies and V is a typical velocity, the relative 
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magnitude of the terms in (4~3) may be estimated by replacing 
V with L 1 and v by V; vector character is neglected in this 
crude approach. It then follows that the right hand term 
is negligible only if the dimensionless ’’magnetic Reynolds 
number" 


R m = ^ o ctLV (4-3) 

is much larger than unity. It is important to note that 

increasing the scale of the flow increases this is why 

conditions of "perfect conduction" occur much more readily on 

astrophy sical and geophysical scales than in the laboratory. 

The discussion in this work will be mainly concerned with 

the limiting case R for it is then that the fluid velocity 

m 

satisfies the criterion (l-2) for the velocity of the lines of 
force. In ordinary hydrodynamics, flows with large Reynolds 
numbers are unstable and tend to break up into turbulent eddies. 
While turbulence may also occur in the flow of conducting fluids 
and of plasmas 9 such behavior will not be considered here since 
relatively little has been established about this subject to 
date (see Batchelor s 1950). 

The most important properties of B following eq . (l-2) are 
flux preservation and line preservation , 

A flow v is said to be f lux preserving for a vector field 
Q if the flux of Q through a surface defined by a closed filament 
of fluid particles is conserved. It will now be shown that for a 
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given magnetic field B a flow v is flux preserving if (l-2) 
is satisfied. 

Let a fluid filament initially following the closed contour 
S be given and let $ be the initial flux of B through it. A 
short interval dt later, each element d [1 of the contour will 
have been displaced by an amount v dt, sweeping in the process 
an area (v x 6i) dt (Figure l). In this time interval, $ changes 
by an amount d$ , ascribable to two causes. The time variation 
of the field contributes the surface intergral 

fl . dA dt (4-4) 

i st 

while the variation of the area bounded by the filament adds 
the flux through the area swept by it (figure l) , equaling 


j> x (v x dX) dt - — J' 7 x (v x B) x dA dt (4-5) 


Thu s 


d$ = f 


( (v x B)| 


x dA dt 


(4-6) 


and this vanishes if (l-2) holds. Because the contour is 
arbitrary, (l-2) is not only sufficient but also necessary for 
flux preservation. For a nonsolenoidal vector field Q the 
condition for flux preservation becomes 

SO 

r= - V X (v x'a) + v (I • a) = 0 (4-7) 
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A motion is called line preserving if two fluid particles 


initially connected by a line of force continue to be so 

throughout their motion . A flux-preserving motion is always 
line preserving, by the following 

argument (for an alternative derivation, see Newcomb, 1958)* 
Consider a sheet ABCD of fluid particles which forms part 
of a magnetic flux tube (Figure 2), with points A and B chosen 
so as to lie, at some initial instant, on the same line of force. 
The magnetic flux through the sheet and through any closed contour 
on it then vanishes, and by flux preservation this property is 
maintained as time progresses. The sheet will therefore always 
be part of a flux tube. 

The same arguments hold for a sheet ABC 1 D* on another flux 
tube which initially intersects the first one along the field 
line AB . As time passes, A and B will continue to belong to 
both flux tubes and therefore lie on their intersection, which 
defines a line of force. Since the choice of A and B is in no 
other way restricted, the motion is line preserving. 

The condition for line preservation is less strict than that 
for flux preservation. By the Helmhol t z-Zorawski criterion 
(Zorawski, 1900; for further references see Truesdell, 1954°> 
Truesdell and Toupin , i960; also Newcomb, 1958, eq . 85) a nec- 
essary and sufficient condition for line preservation is 

x (v x g) + v (v • 2)jxg=0 (4-3) 

The above demonstrates that flnx preservation is not necessary 
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for line preservation; for instance, motion in a force-free 
field with finite, constant resistivity is line preserving but 
not flux preserving. The following proof is due to Prim and 
Truesdell (1950). Let the equation of an open fluid filament 
at time t, given in terms of some parameter s (e.g. the arc 
distance from a given particle in the filament) be 

x = x (s, t) (4-9) 

Then dx/ds is tangent to the filament and the condition 
for line preservation is that if initially 

M X Q = 0 (4-10) 

ds 


so that the filament is aligned with a field line of Q, then 
this condition is met at all times. The time variation of the 
quantity under consideration is 


£(iH ■ (g *2) qts-s) 


dx \ 

51 ‘ 7 2l) x a + 


?>x am 

9S x dt 


(4-11) 


Now ( 4- 10} is equivalent to 


3x 

S’ xs 


(4-12) 


If (4-8) holds and (4-10) is initially given, one gets for the 
initial variation of the expression on the left of (4-10) 
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A(!| k2 ) =x 2 x (s- 2 .vv) 


(4-13) 


= x | dir “ v x ( x x fi) ^ v (v • <g)j = o 

and therefore (4” 10) holds for later times as wello 
(5) Cauchy ,fl s Integral 

Equation (l~2) for the motion of a perfectly conducting 
fluid involves both B and v. For those cases in which both 
quantities have to be derived theoretically another equation 
involving both is required, as well as a set of appropriate 
boundary conditions . The additional equation is the hydro- 
magnetic generalization of the Navier -Stoke s formula and the 
solution of the combined set, even in simple cases, is a for- 
midable tasko 

In practice, however, v may often be obtained separately, 
as for example, when it may be deduced from observations (e.g. 
the differential rotation of the sun) or when magnetic forces 
are negligible in comparison with others acting on the fluid, in 
which case the hydromagnetic equation may be separately solved 
for Vo In such cases one may regard v in eq a (l-2) as given 
and seek a solution for B e 

As was noted before, eq 0 (l-*2) was first studied with Q 
replacing B s in which case it describes the development of 
vorticity in an ideal fluid. It was in this context that its 
general solution was achieved by Cauchy (l8l5)> who transformed 
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the flux-conservation of Q into the material conservation of a 
related quantity. 

Equation (1-2) may be rewritten 

_ |i + v . vb _ B • Vv - B (V • v) (5-1) 

at o l — — — 

Substituting the equation of continuity 



have to be introduced. Its derivation below follows Elsasser 
(1956); some other derivations may be found in the works of 
Brand (1947; §123), Lamb (1879; §146), Goldstein (i960; 

§4.2, 4.3), Truesdel 1 ( 1954) and Truesdell and Toupin (i960). 

Most of these describe the development of fluid vorticity; the 
first application to magnetic fields in perfectly conducting 
fluids is due to Lundquist and Cowling, (Lundquist 1951? 1952) e 
The material description of fluid flow involves referring 
variables to their values at some initial instant t - 0; such 
initial values will be distinguished by superscript (0) c Con- 
sider an arbitrary fluid particle; in an orthogonal cartesian 
reference frame its position vector x at any time t is uniquely 
determined by its initial value x^°^ and by t 

X = X (x^°! t) (5-6) 

The Jacobian J of this transformation between x4^ and x (t being 
assumed given) is the ratio by which a volume element moving 
with the fluid has changed from its initial value, and therefore, 
the ratio between initial and final densities in the particle^s 
vicinity 


J 




Xl , X2 s X 3> 




,(°) „ (°) 


X2 


x 3 


__ <£> 

( 0)1 P 


(5-7) 


Consequently, as long as this ratio is finite, a unique inverse 
transformation exists 

x^°) - x^°^ (x, t) (5“8) 
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The two transformations are related, for any given t, by 


= 6 u 


(in the last equation and in those following the summation 
convention for repeated indices is used) . 

From this, whether i equals 1 or not, 


taj' = Sx k + a^ ) dt U-M 


Writing (5*3) in component form and substituting (5-10) gives 


/ B.\ B i dx 1 ( \ 

(f)+7 47 3 F V^t) = 


It now becomes evident that the integrating factor is 
d^/dx^ Multiplying by it, summing over i and using (5-9) 


leads to 


d / B i ^i } 


dt p ax. 


and since 


Sx . t-o 
1 



I 


(5-12) integrates to 


(o) D (o) 6^j 

B i = & B. _ 

3 p 1 ax. 


(5-13) 


Usually, B^ is given and B^^ required; using eq„ (5-9) then 
gives Cauchy's celebrated integral (Cauchy, 1815) 


B. = p R (o) 1 

1 ^ j a? 


(5-U) 


To investigate the form of eq. (5-14) in curvilinear coordinates 
we adopt tensor notion and indicate contravariant components 
of vectors by upper indices. 

Equation (5-14) > in cartesian coordinates, is then written 


B 1 fL R (°)j m! 

’ p ( °> 


(5-15) 


Let this be transformed into curvilinear coordinates y 1 . 
The transformation of is obtained from the chain 

ru le of differentiation 


_dx 


= dx 


dx 


(o) j 


Jr 


54 


(o) 1 


dy dy 


(°) 1 Sx (o) j 


(5-16) 

| 


while B 1 and B ° ^ transform as contravariant vector components., 
Substituting into (5-15) and using (5-9) > which is valid for the 
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y as well, yields an equation having the same form as (5-15). 
That equation is therefore valid for all x 1 , rectilinear and 
curvilinear. In particular, in orthogonal curvilinear 
coordinates with scale factors h^ (e.g. Morse & Feshbach, 1953> 
§1.3) the proper components of B (i.e. those defined in 

terms of unit vectors) satisfy 


B (i) = 


_P_ h. 

» 1 


(5-17) 


(summation on j only) 

Cauchy's integral is frequently written 


_P_ fi (o) . v (o) x 


( 5 - 18 ) 


with similar notation in the equations leading to it. The 
invariance of eq. (5-15) under coordinate transformation shows, 
however, that V^°^x does not transform in the manner of a 
second order tensor such as Vv [Such a transformation was 
assumed by Elsasser (1956), leading to the incorrect eq. (7.12) 
in his work]. Strictly speaking, V^°^x is a double tensor field 
(Michal 1927, 1947; Truesdell and Toupin, i960; Ericksen, i960), 
the components of which may depend on two points in space - here 
x and - and the transformation of which may simultaneously 


involve the local vector bases at both points. In the trans- 
formation of 7^°^x, for instance, transforms at x^°^ while 


x transforms at the particle* s given location; this is expressed 
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by eq„ (5-l6) o For a more rigorous and comprehensive treat- 
ment of double tensor fields and their applications the reader 
is referred to the works of Truesdell and Toupin (i960) and of 
Ericksen (l960) o 

Cauchy’s integral may be generalized for finite conductivity 
(Mauarsberger , 1964) ° Let E* be the electric field in a frame 
moving with the fluid, which to first order in v/c equals 

i 

Ey — E + ' (v x B) 


Then ( 4“ 1 ) becomes 


M -V x (v X B) = - ~(V X E*) 

dt “ 


(5-19) 


Retracing previous calculations, one obtains 


A. • Vv = - - (v X E* ' 

dt V p / p ~ P 


(5-20) 


and with the same integrating factor as before this gives 


_d_ / El • Vx^ 9 2 \ = - 1 
dt \ P / P 


1 (V x E*) • Vx 1 

P “ 


(o) 


( 5 - 21 ) 


Integration then results in 


B(x) = 


£_ f b ( °^ 
(o) k 


J fc j(V x E*) • v4°4t 


i- V (o) x (5-22) 
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where the integration extends along the track of that particle 
which has position vector x at time t. The Jacobian J is 
defined by eq. (5-7) and may be expressed in terms V^°^x 
(compare Truesdell 1954, eq. 84-l) , which is done in Maiersberger ' s 
work. Mauarsberger also derives the magnetic analog of Ertel's 
vorticity formula and some applications; for more details, the 
reader is referred to his article. 

(6) Application of Cauchy’s Integral 

Let dx be the (time dependent) eLement of length occupied by 
an infinitesimal fluid filament. If at t = 0 the filament is 
aligned with the magnetic field, one may write 

dJ°) = B^dX (6-1) 

with dX an infinitesimal constant. Using the chain rule of 
differentiation on (5-6) and the notation of (5-18) gives, for 
t > 0 

(o) 

dx - dx^°^ . V^°^x - —— — B dX (6-2) 

The above implies that the motion is line -preserving as could / 
of course^ be directly inferred from (4-8). Conversely, Cauchy’s 
integral may be easily derived from the line-preservation and 
f lux-preservation properties of B which follow from (l-2) 
(Lundquist, 1951, 1952; result is expressed there in terms of 



Eliminatind dX between the absolute values of (6-1) and 
(6-2) yields 
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(6-3) 


B 

B ( °) 




from which the intuitive interpretation of Cauchy's integral 
may be summed up in three statements: 

(1) The motion is line pre serving . 

(2) The magnetic field intensity is amplified in direct 
proportion to the stretching dx/dx^ of the fluid 
along the field's direction. 

(3) The magnetic field is amplified in proportion to the 
bulk compression p/p^°^ 0 

Most qualitative discussions of phenomena involving motion 
of lines of force are based on one or more of the foregoing. 
Consider the differential rotation of the sun: if initially 
the lines of force of the solar magnetic field lie in meridional 
planes, then the differential rotation of the solar surface, 
because of line preservation, would tend to wind them into a 
pair of spirals having opposite sense in opposite hemispheres. 
The same process stretches solar matter along lines of force 
and may therefore be expected to amplify the field. This is 
roughly the basis of Babcock's theory of the solar cycle (l96l); 
the energy required for the field's amplification is ultimately 
supplied by the mechanism which maintains the differential 
rotation . 

Again, when a star condenses out of a large volume of 
tenuous conducting gas permeated by a weak magnetic field, the 
field is likely to undergo considerable amplification due to 
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compression. Since the field's final strength maybe high enough to 
contribute appreciably to the forces acting on the gas, the 
process is important in theories of stellar evolution. (Hoyle, 

1957 ) . 

Equations ( 6 - 2 ) and ( 6 - 3 ) may also be used numerically. By 
following a small fluid element from a point at which B is known 
and noting the variation of its density and its dimension along 
the field, B at other points along the element's motion can be 
found . 

Finally, situations exist in which eq. (5-18) may be used 
analytically. Consider, for instance, the solar wind emanating 
from the rotating sun. It can be assumed (Stern 1964 ) that near 
the sun the lines of force are radial (if one starts from a 
rigidly rotating sun problems of continuity arise) . Assuming 
that the flow velocity v is radial and that both it and the 
angular velocity to of the field's source are constant, the con- 
travariant components of eq. (5-6), in spherical coordinates 
corotating with the sun, are 


(o) 

r = r 4* vt 


e = 

A°) nt >) 

0=0 — Ut = 0 


0 ) / (oh 
— ( r — r ) 
v ' ' 


( 6 - 4 ) 


Assuming also 



( 6 - 5 ) 
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equation (5-17) gives 



leading to the well-known angle (Parker, 1958) 

X=arctg (u) r sin 0/v) (6-7) 

between the radial direction and that of B. Another example 
of this type, an application of Cauchy 1 s integral to a problem 
in dynamo theory, has been given by Parker (1955)° 

(7) Neutral Points in the Magnetic Field 

We now consider the case in which the conductivi ty a is finite, 
though still very large c When the velocity v is negligible, eq» 
(4-2) reduces to 

M = — v 2 b (7-1) 

St M- 0 cr - 

which yields a parabolic diffusion equation for each of the 
cartesian components of B 0 The field then tends to decay ex- 
ponentially with a decay time of order T — p, aL 2 , where L is a 
typical dimension of the field; the decay of the flux passing a 


given volume can also be visualized as a diffusive spreading-out 
of the lines of force with velocity of order L/t =(jj, Q aL) 1 . An 
example of this is the decay of the magnetic field due to a 
current distribution in a conducting sphere (Elsasser, 1956). 

The foregoing suggests that the general solution of (4-2) 
is a blend of convection of the field by the fluid's velocity v 
and its diffusive decay, the relative importance of the two effects 
being determined by R . A useful concept in this case, introduced 
by Sweet (1950), is the slippage velocity of lines of force. 

Sweet visualized the velocity of lines of force in a conducting 
fluid as consisting of two components: a convection with the 

fluid's velocity v and a slippage w due to finite conductivity. 

In addition to finite conductivity he also considered the terms 
added to E# in Ohm’s law for dense plasmas (_Spitzer, 1962 eq. 

2-21; the largest of these terms appear in eq. (1-5)] lumping 
them together in an ''impressed e.m.f." Since the total 

velocity of the lines of force has to satisfy (l-2), one obtains 

Vx (wxB)+E = V 2 B (7-2) 

~ext gL Q cr “ 

The component of w parallel to the field is undefined in (7-2) 
and it is natural to add 

(w ■ B) = 0 (7-3) 

In what follows E. ext will be neglected. Removing the curl and 
isolating w then gives 
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(7-4) 


where v/ 1 is a relabeling velocity satisfying eq, (3-4) . 

Disregarding wj , the slippage velocity is proportional to the 
magnetic body force acting on the fluid. One thus gets the 
picture of a viscous interaction between the fluid and the 
lines of force, with the lines exerting a force on the fluid 
proportional to their relative velocity. The order of magnitude 
of w is 1 and it thus furnishes a more precise definition 

of the "diffusion velocity" previously introduced. 

The analogy with viscosity can be expressed in a more 
quantitative manner. It may be shown (for details and applications, 
see Batchelor, 1950, and Landau & Lifshitz, i960, §55) that in an 
incompressible fluid with viscosity y, subject to gravity and 
pressure, equation (l-l) for the vorticity Q is replaced by 


^ “ V x (v x q) — i;V 2 q 

which formally resembles (4-2), Because of this analogy, the 
qu an tity 

U m = W' 

is sometimes called the magnetic viscosity (Elsasser, 1956, 
eq. 2,13) , 

Equation (4-2) is o^ prime importance in the theory of 
self-sustaining fluid dynamos. Special cases have been solved 
by Elsasser (1946a, 1946b, 1947), Parker and Krook (1956), 
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Stern (1964) and others, but the general treatment is difficult; 
it does not reduce to scalar diffusion equations in the convected 
frame for a and B • This is in keeping with the general rule 
that, except for simple cases, equations involving the vector 
Laplacian do not reduce to a set of corresponding equations 
involving scalar Laplacians. 

In astrophy sical applications many cases occur in which <7 
is exceedingly large. It has been pointed out by Sweet (1956) 
and later by Parker (1957) and Pestschek (1963) that in these 
cases, if the magnetic field possesses neutral points at which 
its intensity vanishes, dissipative effects can be greatly 
accelerated in the vicinity of these points. Following these 
authors, we now investigate highly conductive fluids in which 
the effects of finite conductivity are small everywhere except 
in the neighborhood of neutral points. 

Astrophy sical interest in neutral points dates from an 
observation by Giovanelli (1947) that solar flares occur pre- 
ferentially in the vicinity of neutral points of sunspot fields. 
This idea - which still lacks complete confirmation because of 
the difficulty in identifying neutral points on the sun - led 
Dungey (l953> 1958, 1963) to develop a theory of acceleration of 
charged particles by electric fields near neutral points. Much 
of what is said here about the classification and behavior of 
neutral points follows Dungey 1 s work. 

The magnetic field B near a neutral point N may be expressed, 
to lowest order, by a Taylor expansion 
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( 7 - 5 ) 


B = a o r 

where r is the radius vector to N and a " 1 = VB|^ (superscript T 
denotes transpose) is a constant dyadic with zero trace, 
symmetric if the current density at N vanishes ■ The direction 
of B at N is undefined; nevertheless, it is useful to define 
a n line of force through N" as a continuous smooth line passing 
N and tangent to B at all other points (the term "neutral line" 
is sometimes used, but this is apt to be confused with a line 
along which B vanishes) • 

Neutral points can be classified (Dungey 1953) by the number 
of lines of force passing them. Let such a line have a unit 
tangent £ at N, At points on the line infinitesimally close 
to N 

r=)i[ 

B = \a - £ = |i(\) f. 

the last equality following from the condition that B is tangent 
to the line. Since X and p, are both real (though infinitesimal), 
£ must be an eigenvector of a corresponding to a real eigenvalue. 
A line though N on which B vanishes is also considered here a 
line of force, corresponding to an eigenvalue zero. A bundle 
of lines of force having a common tangent at N is counted as one 
line only. 

There exist therefore two kinds of neutral points. If a has 
only one real eigenvalue, there is only one line of force passing 
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through N, which is then called (by the shape of lines of 
force in its vicinity) an O-type neutral point. The behavior 
of B near an O-type neutral point forms the basis of Cowling's 
theorem (Cowling, 1934 ) on the impossibility of stationary 
axisymmetric fluid dynamos. Alternatively, there may exist 3 
real eigenvalues X , with corresponding unit vectors £ , in 
which case 3 lines of force or bundles of lines pass through 
N and it is known as an X-type neutral point. All current-free 
neutral points belong to this class, for a is then symmetric 
and has only real eigenvalues. In what follows, only X-type 
neutral points will be considered. 

Because the divergence of B vanishes, the three eigenvalues 
add up to zero and cannot all be of the same sign. It will be 
assumed here that Xi and X2 are positive and X3 negative; the 
alternate case of two negative eigenvalues is then obtained by 
reversing the magnetic field everywhere, a procedure which, it 
will be noted, does not change the disposition of the magnetic 
force. The vectors £^ will be chosen of unit length and their 
sense is determined by stipulating that the scalar products 
(Ii * &)> (£2 • 13) and (lila S 3 ) are positive [the sign of 
(£1 ' £3) cannot be freely chosen]. 

If r lies in the plane of £ 1 and £2 » B is contained in that 
plane as well. Choosing an orthogonal cartesian system for which 
this is the (x,y) plane, it is found that everywhere in this plane 
the quadratic form 

r • B = a:r r 


- 32 - 



involves only the top left (2x2) submatrix a f of a, which has 
the positive eigenvalues X^ and X 2 , If in addition no current 
flows in the z-direction, a" is symmetric as well and therefore 
positive definite. Then ^ and 4,2 are orthogonal, r * B is 
positive for all directions and the lines of force point outwards 
everywhere in the (£1 £ 2 ) plane, In this case each of the lines 
of force to which Ej and £2 are tangent spreads out into a bundle 
of outward-going lines, the two bundles, between them, covering 
the plane. In the general case in which £1 and £2 are inclinced 
at some angle 0, a^ is positive definite only if 

cos 2 e < 4 X 1 X 2 / (X! + \z) Z (7-6) 

In other cases the lines of force diverging fron N bend so as to 
approach N again before finally heading away (compare Dungey, 
1963) . 

The claim is often made that lines of force are broken and 
reconnected at neutral points. This will now be examined. 

In the approximation (7-5) the cartesian components of B are 
linear combinations of those of r and therefore, in this order 
of approximation, V 2 B vanishes. Neglecting E ^ in (7 = 2), the 
slippage velocity w then satisfies (3-4) and therefore has the 
character of a "relabeling velocity" , by which lines of force in 
a given pattern are transformed among themselves, without change 
to the pattern as a whole. 

Assume for simplicity a field in which there is no fluid 
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motion: in such a field, w will be the entire velocity of 

field lines. Near the neutral point, w is large (it diverges 
at N) and it transforms field lines within a fixed pattern, but 
at large distances fron N it diminishes rapidly to a negligible 
magnitude- 

Consider a line of force which initially passes points ABC 
near N (Figure 3). A time interval dt later the line will have 
shifted a distance w dt. At the "roots" of the line far from N 
this shift is negligible, so that if one defines a line by its 
roots one will be still observing essentially the same line. 

Near N, however, there will be a real change, with the points 
ABC moving to DEF and the line becoming differently connected. 

The slippage velocity becomes infinite at N. This is necessary 
to enable point A to "turn a corner" at N and transform to D 
while keeping up with B and C as they transform smoothly into 
E and F. As a consequence, no matter what the flow velocity v 
of the conducting fluid is, there will always exist a small 
neighborhood of N in which the convection of field lines by the 
fluid velocity can be neglected and their velocity be taken as w. 

To obtain some more insight into the slippage process near 
neutral points we introduce an essentially two-dimensional model, 
with vectors Ej and £2 making an angle 0 in the (x, y) plane and 
£3 following the z-axis, and with corresponding eigenvalues X, — X 
and 0. The diffusion velocity w along the direction of (^ 4 - £ 2 ) 
is then radial (this corresponds to diffusion fron N to D in 
Figure 3) and its magnitude at distance x is 
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cos 0 


( 7 - 7 ) 


I ^ ~ ~\^~o x sin 1(0/2) 

One notes that the magnitude of B s represented by \ 9 is 
absent and that w increases as the angle between the axes 
diminishes (it is zero for 0 s ft/2 because then V x B vanishes) 
and also is inversely proportional to the distance from N» 

(8) The Evolution of Neutral Points 

The volume force on the fluid (T denotes transpose) is 

F 1 j x B * — (a — a^ ) 0 a 0 r (8-l) 

“ l^o 

The force vanishes at the neutral point and, because of its 
orthogonality to B, has no radial component on the main axes 
defined by the ^ « In general , the radial component of F at a 
point 

r * r u. 
is given by 

r F r - | 2 (SjSj) u. Uj (x. ~\.) 2 (8-2) 

This component vanishes on the quadric surface obtained by 

equating (8-2) to zero* which is evidently a cone* since it 
contains three non-coplanar straight lines (the main axes) 
intersecting at the origin Q The force then pushes inward 
inside the cone and away fron N outside it -■ or vice versa, 
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depending on the sign of • £ 3 ). If the fluid flow follows 

this force, it will tend either to open the cone up into the 
vicinity of a plane ("the neutral plane 11 ) or to close it into a 
line (’’neutral filament") , either of them passing through N. 

Such "folding” of the field was obtained by Dungey (1953, 
1958 ) when he first treated the fluid's motion near the neutral 
point* His calculation assumed that pressure and conductivity 
can be neglected and that the fluid velocity vanishes at the 
neutral point, so that it can be expanded, analogously to ( 7 - 5 ) 

v = b • r ( 8 - 3 ) 

Equation (l-2) then gives a (dyadic) equation 


-= = (b • a) —(a . b) — a Eb . . (8-4) 

o t 1 1 

(forming the scalar product with r recovers 1-2), Similarly, 
the force equation 

!“ + v • Vv = (j x B) (8-5) 

reduces to 

+(b • b) = — (a — a T ) • a ( 8 - 6 ) 

at 

If initially both v and j vanish, so that the magnetic axes are 
orthogonal and a is symmetric, Dungey showed that a small perturbation 
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amplifies itself and produces a flow diminishing the angle 

between axes of oppositely signed eigenvalues (angle 0 in 

Figure 3)- The current density j meanwhile increases, until 

the electric field E* (in the fluid's frame) becomes appreciable. 

Dungey called this an M instability ,, ; it was his idea that such a 

development takes place in solar flares and that the current 

then becomes a "discharge" accelerating particles to high energies. 

While Dungey ! s derivation is formally correct, caution must be 
exercised in its interpretation, for the growth of a local instabil- 
ity associated with N contradicts the existence of stable neutral 
points*. Such points do exist*. Consider for instance two fixed 
parallel magnetic dipoles, directed along the line connecting 
them and immersed in an infinite homogeneous conducting fluid at 
rest (Figure 4)- In such a conf iguration a neutral point (more 
precisely, a circular locus of neutral points) is formed between 
the dipoles, resembling in all respects the one investigated by 
Dungey «* If this system were unstable, one would still expect that 
after all fluid motions and volume currents had decayed there would 
remain a homogeneous medium with two embedded dipoles. The final 
energy then equals the initial one and there is none left to drive 
the instability. 

To understand Dungey 1 s results one notes that the effect deduced 
by him is not confined to the neighborhood of N. Indeed, the 
folding of the main axes and of the associated system of field 
lines extends as far as the approximation (7-5) does, and it is 
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not at all certain that some surface can be found, enclosing 
a relatively small region around N, beyond which the perturbation 
is negligible. In the example given, at least, it seems that 
no such surface exists and that the "folding 11 near N is only 
part of a deformation involving the entire field, because the 
only source of energy available here is the relative potential 
energy of the two dipoles. In particular, allowing the dipoles 
to attract each other (Jackson 1962, §4.2) through a finite 
distance will release energy and will also fold the lines near 
the neutral point in the manner described by Dungey. It thus 
seems that Dungey 1 s mechanism is actually the description of a 
lack of equilibrium in the entire field rather than of an 
instability associated with the neutral points. 

A problem somewhat similar to Dungey* s was treated by Chapman 
and Kendall (1963) who considered a line along which the field 
vanishes, surrounded by a tube of conducting fluid with finite 
cross section • In this case Dungey* s instability proceeds without 
reference to the magnetic sources, since the fluid is decoupled 
from them by the vacuum surrounding the tube. Such decoupling 
seems unlikely, however, in astrophy sical applications. 

To obtain any local behavior at the neutral point while the 
field is held fixed at some distance from it, one has to extend 
the expansion (7-5) at least to second order. In principle, 
equations (8-4) and (8-6), generalized to this order, or the use 
of an energy principle, would then show under what conditions, if 
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any, can a purely local instability proceed at the neutral point = 

Attention was first called to the importance of motion of the 
field's sources by Sweet (l956 3 1958) „ He considered in particular 
the neutral point between the two dipoles described previously 
(Figure 4) and assumed each dipole to be embedded in a mass of 
conducting fluid, the two masses moving in opposite directions 
towards the neutral point „ 

Sweet did not restrict the model to electromagnetic forces 
alone, In this connection, he noted that forces due to pressure 
gradients will predominate near the neutral point, since the 
electromagnetic force tends there to zero. The flow produced 
by mechanical forces when two fluid masses are pressed together 
is similar to the one occurring in Dungey 1 s model: the fluid is 

squeezed out along a plane or line. If the neutral point is 
initially on such a "neutral plane" (or possibly, "neutral filamen 
Sweet expected the flow to distort the magnetic field as was pre 
viously described, leading to its flattening near the plane and 
to the build-op of current density and electromagnetic forces, Th 
process is finally limited by the necessity of maintaining a narro 
region by which the fluid may escape from the neighborhood of the 
neutral point,,. It is in this final state that enhanced diffusion 
near the neutral point takes place. 

The following simple model is due to Parker, who followed Sweet 
ideas quan titatively in considerable detail (l957 s 1963) « Parker 
assumed that the region around the neutral point N collapses into 
a thin "boundary layer" of constant thickness 26. The fluid 


approaches the layer with velocity u along the direction 

normal to it, which will be chosen as the x axis; the y direction, 

along which the fluid is squeezed out, is then that of the field 


far from the layer, which is assumed to be of constant 
magnitude but having opposite polarities on different sides of 
the layer. The magnetic pressure on the flanks of the boundary 
layer is balanced, near N, by an increase in fluid pressure. 

Let the layer have a finite length 2L in the y direction. At 
| y| = L, the fluid leaves the layer and if the field and pressure 
have dropped close to zero near this point, the fluid's kinetic 
energy density p v 2 /2 will be roughly that acquired by a pressure 
drop of B q 2 /2|jl o * This means, then, that the velocity of the fluid 
leaving the layer is of the order of the Alfven velocity 




(8-7) 


By continuity 

V A 6 = W. L (8-8) 

and if u and L are given, 6 may be found. 

Inside the layer opposing field lines diffuse towards each 
other with a velocity of the order 6) 1 and ultimately merge. 

Outside the boundary, field lines move with the fluid's velocity 
U. In order for the two motions to be continuous, one needs 
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u 


( 8 - 9 ) 


« 6)' 1 

Eliminating 6 between the last two equations gives 


u « (V^o a L) ^ (8-10) 

In applying the above to the sun Parker assumed that a solar 
flare arises when two regions containing opposing magnetic 
fields are pushed against each other with velocity u e The 
field intensity; density and the total energy released can be 
estimated from observations; from the total energy* the thickness 
D of the region of field annihilation is found to be of the order 
of 10* km* The total time D/y k for the annihilation of magnetic 
energy is then 10* - 10 6 seconds, as compared to observed times 
of order 10 seconds 0 

The above argument may be modified taking into account more 
details of the field c s structure near the neutral point* assuming 
that it has the form given in Figure 3 with the angle 0 small but 
finite c. There will exist a region around N in which w predominates, 
and as a simple approximation one may neglect v in this region and 
neglect w outside it. If 6 is the extent of the region in x 
direction (along the line NA in Figure 3) then its length in the y 
direction will be of the order 6/sin ^ 9/ 

Equation (8-9) then expresses the fact that on the boundary 
of this central region, on the x-axis, the slippage velocity equals 
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the impinging fluid velocity. However, since this transition 
occurs close to a stagnation point of the flow, the impinging 
velocity may be much less than the velocity u of the fluid far 
from the layer. Writing the impinging velocity as eu , where 
e < 1, one finds that a factor e ^ must be added to equation 
(8-10). If e is very small, this increases u considerably and 
achieves a corresponding reduction in annihilation time. 

In addition, unless the angle 9 is extremely small, the 
central region will occupy only part of the boundary layer, 
terminating at about | y| = 6/sin £ 9- Beyond this point the 
layer is likely to broaden; it will certainly do so if the flow 
resembles that given by equation ( 8 - 3 ). Equation ( 8 - 8 ) then also 
requires modification and again this leads to an increase in u. 

A more elaborate extension of Parker 1 s model was postulated 
by Petscheck (1963), who divided the boundary layer into two 
regions. In a section 0 < | y| < y* straddling the neutral point 
the mechanism suggested by Sweet and Parker holds; for y*<| yj < L, 
however, the layer diverges at an angle of about 0.1 radians, the 
fluid velocity is maintained near and magnetic energy is 
removed by propagation of hydromagnetic waves. In application to 
solar flares, times of the order of 100 seconds are obtained, 
because the length 2 y# of the central region is now much less than 
the total length 2L and therefore the thickness 20 near N is 
corresponding ly reduced. 

Petschek also derives an upper limit to the rate of line re- 
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combination, The main disadvantages of his approach are 
that it is not clear that the "wave region" envisioned will 
actually be created and that his assumption of a curl -free 
magnetic field everywhere outside the boundary layer leads to 
the correct maximum recombination rate„ 

In cone lii sion it should be stressed that the theory of 
neutral points is very far from complete at this stage Q Almost 
all existing theories assume two-dimensional structure, the 
achievement of a nearly time -independent collapsed state and 
the stagnation of the fluid flow at the neutral point, any or 
all of which may not be justified,, Even if an equilibrium in 
the collapsed state occurs, it may be unstable (Parker, 1963; 

Jaggi, 1963; Severny, 1963) • Whether such instability merely 
leads to small-scale turbulence and to an increase in effective 
resistivity or actually disrupts the structure of the boundary 
layer remains yet to be seen 0 

The most important deficiency of the theory of neutral 
points, however, is the lack of experimental conf irmation , 

Though the theory was originally conceived to explain solar flares, 
the complexity of sunspot fields and of flares and the difficulty 
in observing them have prevented, so far, any conclusive proof 
of an association between flares and neutral points. There may 
be some additional evidence from the "tail" of the earth's magneto- 
sphere (Ness, 1965) where a region of weak magnetic field, commonly 
called "neutral plane", has been observed near the magnetic equator 
on the night side of the earth, separating oppositely directed 
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fields. Both the magnetic field and particle density are much 
lower here than near solar flares; whether rapid reconnection 
of field lines (Axford et al.,1965) and particle acceleration 
(Speiser, 1965; Coppi et al., 1966) occur under these circumstance 
is still subject to controversy (Dessler, 1965). 

Finally, an attempt to investigate neutral points experi- 
mentally in the laboratory has been recently started (Bratenahl 
& Hirsch, 1966). The neutral point (or rather, a line of neutral 
points) was obtained by allowing two expanding cylindrical fields 
to collide, leading to a situation somewhat resembling that in- 
vestigated theoretically by Green (Green, 1965). Observations 
of the evolution of these fields strongly suggest that a recon- 
nection of magnetic field lines does in fact occur. However, 
instead of the field tending to a relatively slowly varying 
collapsed state, oscillations are observed with field lines near 
the neutral point alternately folding and opening. 

It thus appears that what is most needed in the study of 
neutral points is observational evidence, either in the laboratory 
or in space. Until such evidence becomes available, the preceding 
can serve as no more than a tentative theory. 
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- No caption 

- Two Interesecting Flux Tubes 

- The Magnetic Field Lines Near an X-type Net'tral Point 

- The Neutral Point Formed Between Two Parallel Dipoles 
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